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^' Abstract 

An analytical approach is carried out that provides an inviscid stability criterion for the 

(N ■ 

strato-rotational instability (in short SRI) occurring in a Taylor-Couette system. The control 

<n : 

parameters of the problem are the rotation ratio ii and the radius ratio rj. The study is 
motivated by recent experimental [1] and numerical [21 [3] results reporting the existence of 
unstable modes beyond the Rayleigh line for centrifugal instability (fj, = rj 2 ). The modified 

o 



Rayleigh criterion for stably stratified flows provides the instability condition, \i < 1, while 
in experiments unstable modes were never found beyond the line ji = rj. Taking into account 

•i-H , 

finite gap effects, we consider non axisymmetric perturbations with azimuthal wavenumber 
I in the limit IFr << 1, where Fr is the Froude number. We derive a necessary condition 
for instability : jjl < fi* where //*, a function of rj, takes the asymptotic values, fj,* — ► 1 in 
the narrow gap limit, and jx* — ► 2tj 2 {1 + rj) in the wide gap limit, in agreement with recent 
numerical findings. A stronger condition, \i < rj, is found when rj > 0.38, in agreement with 
experimental results obtained for rj = 0.8. Whatever the gap size, instability is predicted for 
values of ji larger than the critical one, ji c = r/ 2 , corresponding to centrifugal instability. 



1 Introduction 

Despite its long-lasting study since the work of Taylor [4] , the stability of cylindrical Couette 
flow remains a vivid research area attracting many investigators. A survey of the literature 
on the topic can be found in [5j. The flow occurs in the annular gap between two concen- 
tric cylinders of radii R\ < R2 rotating independently at the angular velocities 0,± and O2, 
respectively. The control parameters are the radius ratio r\ = R1/R2 and the rotation ratio 
\i = f^/f^i. Inside the gap, the angular velocity profile of the laminar steady flow is O(r). Its 
linear stability with respect to axisymmetric perturbations is governed, in the inviscid limit, 
by the Rayleigh criterion [6] : d(r 2 VL) 2 / dr > 0. A generalized Rayleigh criterion for non- 
axisymmetric centrifugal instabilities has been derived for a free axisymmetric vortex using a 
large axial wavenumber WKB approximation [7J. The instability takes the form of a spatially 
oscillating mode localized between two turning points where it matches with the exponentially 
decaying solutions outside. The generalization to take into account a background rotation 
and a stable stratification did not reveal fundamental changes in the Rayleigh criterion. This 
contrasts with what was found for bounded flows [U [10]. In these studies, it was shown 
that a stable stratification drastically changes the stability criterion, showing a strong analogy 
with the modified Rayleigh criterion obtained when a magnetic field parallel to the cylinders' 
axis is present. In the hydromagnetic case [111 [T2] the condition for stability with respect 
to axisymmetric perturbations is : d(Q 2 )/dr > 0. The presence of a stable stratification has 
a similar effect on the Rayleigh criterion though its validity is in that case restricted to non 
axisymmetric perturbations [8l l9l HO]. 

In the astrophysical context, the magneto-rotational instability (MRI) was recognized as a 
potential source of turbulence in accretion disks [13], especially Keplerian disks, with angular 
velocity O(r) ~ r ~ 3 / 2 5 unable to sustain centrifugal instability according to the standard 
Rayleigh criterion. The MRI mechanism is robust, it occurs in bounded as well as unbounded 
flows and for compressible or incompressible fluid. The capacity of the strato-rotational 
instability (SRI) to trigger turbulence in Keplerian disks raises the open question of what 
are the appropriate radial boundary conditions at the edges of the disk. The existence of 
SRI has been demonstrated in the inviscid limit for perturbations that satisfy no normal 
flow conditions on the channel walls. In the viscous case, both the no-slip [2] and stress-free 
|9j conditions give rise to instability. Mixed boundary conditions in which no normal flow 



is imposed on one side of the channel and zero pressure on the other side were not able to 
sustain growing modes [10] . The effect of a vertically varying stratification was also examined 
|10j showing that SRI persists in that case. 

Beside possible applications in geophysics and astrophysics, SRI was studied in the labo- 
ratory. Until very recently the conclusion of experimental studies was that a stable vertical 
stratification stabilizes the flow |14 |, H5| . rT6] . The experimental evidence of the strato-rotational 
instability (SRI) in a Taylor-Couette system was definitely assessed in [Tj. Moreover, the val- 
ues of the control parameters for which instability occurs were found in good agreement with 
numerical predictions |2j that yields the condition : jjl < 77, for instability. 

The aim of the present contribution is to derive an instability criterion by an entirely ana- 
lytical analysis that improves previous studies achieved in the small gap limit [SJ O [TU] . When 
curvature effects are neglected the problem under consideration reduces to the stability of a 
stratified plane Couette flow rotating at constant angular velocity Q. In cartesian coordinates 
(x, y, z) the flow velocity is V = (0, Sx, 0) where S is the constant shear. The condition for 
instability in the stratified case is S/20, < 0, while the standard Rayleigh-Pedley criterion 
gives : 5/2S7 < —1, when there is no stratification. The substitution S — ► rQ , is often 
used to deduce the modified Rayleigh criterion for stratified flows with curved streamlines 
: d(£l 2 )/dr < 0. When applied to the circular Couette flow, the instability criterion simply 
yields : /j < 1. However, in experiments [1] achieved for a finite value of the gap, unstable 
modes were never found beyond the line /j, = rj, in agreement with previous numerical results 
|2j. Quite recently, the situation has changed since numerical calculations [3] have shown 
that the stability line has a more complicated dependence on?). In the narrow gap limit, the 
stability limit was found beyond the line /x = 77, while in the wide gap limit it was found in 
between the lines \i = rj 2 and \x = 77. The narrowing of the instability range, when the gap 
size increases, is possibly due to curvature effects that were not taken into account with suffi- 
cient accuracy in the transposition S — > r$7 . In the present contribution a more appropriate 
treatment of curvature effects is carried out leading to an instability criterion that involves 
the two parameters [i and 77. Our stability results will be compared to the experimental ones 
[Ij obtained for a value of the radius ratio 77 = 0.8, and to the numerical ones [21 [3] for three 
values of the gap in the range 0.3 < 77 < 0.78. 



2 Stratified circular Couette flow 

In cylindrical coordinates (r, <p, z) the velocity field in the basic state is V = (0, r£l(r), 0) with 
the angular velocity given by [2] 



4 + b] = ni«(r) where A = R 2 - — %, and B = ^ — % 
r l J 1 — r] z 1 — r/ z 



The fluid is assumed incompressible with a stable density stratification along the vertical 
cylinders' axis (dp/dz < 0). 

2.1 Linearized equations for perturbations 

The governing linearized equations for the perturbed velocity u = (u,v,w), the pressure and 
entropy perturbations, respectively p and h are given in [Sj. In the study of stratified plane 
Couette flow [H [10] it was particularly convenient to reduce the full set of equations to a set 
of coupled equations for the radial velocity component u and the pressure p. Following the 
same procedure, the perturbed quantities are sought in the form (u,p) = expi(mz + l<p -+- 
ut)[(u(r) , p(r)] where m and I are respectively the axial and azimuthal wavenumbers, these 
notations are consistent with [8] but not with [2] where m has a different meaning. The radial 
dependencies u(r) and p(r) satisfy 

I G 

aD if u + -Zu = —i-prP (2) 

r K 

(a 2 - 2nZ)u - 2i-VLp = -iaDp, (3) 

r 

where D = d/dr and D* = D + 1/r. Here, a = uj — Ifl, and Z = 20, + rO, . The quantities K 

and G are given below 

;2 
K = a 2 -N 2 , and G = a 2 m 2 + — K (4) 

where N is the Brunt- Vaisala frequency. For circular Couette flow with angular velocity 
([1]), the value of Z = 2£liB, is a constant. Thus, the Rayleigh discriminant, <I> = 2QZ, is 
proportional to and it can be written : 

$ = AnjBn = nl$ (5) 

The system ([2])- ([3]) can be reduced to a second order differential equation for the radial 
component u, as in [8]. However, such a formulation leads to tractable results only in a 



few cases, for instance in the large axial wavenumber limit (m >> 1) or the small gap limit 
(77 — > 1). The present stability analysis will closely follow the procedure done in [10] for the 
plane Couette flow. The main difference is the consideration of curvature effects. However, 
the asymptotic analysis |10| based on the small azimuthal wavenumber assumption cannot 
be transposed directly in the circular case where I takes integer values : I = 1, 2, • • •. In the 
circular case, we shall assume the following scalings 

u = in u, Q = Q n, => a = m a, K = -N 2 (l - l 2 Fr 2 a 2 ) (6) 

where Oo = (fii + Q2V2, is the mean angular velocity and Fr = Qq/N, the Froude number. 
The asymptotic value K — ► —N 2 obtained in the limit / << 1 [10] is recovered here in the 
limit IFr « 1, that can be reached either when I « 1 and Fr ~ 0(1) or I ~ 0(1) and 
Fr << 1. Another important assumption concerns the term H = a 2 — 2QZ in Eq. © that 
was approximated in [TO] by H — > — $ where $ the Rayleigh discriminant given in Eq. ([5]) is 
now written $ = ^q$ leading to H = S7q(/ 2 <t 2 — <!). Noticing that # can takes values as large 
as AQ, 2 /Qq, the assumption I « 1 will be replaced by the weaker condition : I < l m ax with 
Imax = 20i/Oo- To illustrate our purpose, this gives the following constraints : I < 2.2 for 
/j = 0.8 and I < 3 for /j = 0.3. Thus, the lower is /j, the higher is the allowed value of I. In the 
following, we shall consider the limit IFr « 1 with I < l ma x- This leads to simplifications 
in Eqs. ([I])-© that are listed below 



K^-N 2 , #->-*, i, = -l 2 [fh 2 * 2 - -)=-l 2 g( r ) (7) 



-• n ^ G=-p( m 2 d 2 -\ 

where m = mFr. 

It must be stressed that the generalization of the Rayleigh criterion [7] performed in the 
limit m >> 1, is not equivalent to the approach described here. The scaling used in [7] for 
the frequency, suitable in describing inertial-gravity waves, could explain why SRI was not 
found in this study. 

2.2 Governing equations in the limit IFr << 1 and I < l max 
The perturbations u and p are scaled according to 

u = lu, and p = 0,Qp (8) 



and substituted in ([2|)-([3]) leading to the following set of equations : 



a D*u H u 

r 

n 

<3?u + 2i— p 



i g p 
i a Dp 



(9) 
(10) 



where £l = Q/CIq and Z = Z/CIq. We have reported in Appendix A the equations that should 
be solved for I > l ma x when the limit H — > — 3> is no longer valid. 

Eqs. ([9j)-(|lUp are the circular version of those obtained for a plane Couette flow [lUj . they 
constitute a set of coupled equations provided <3? 7^ 0. It was shown in [TO] that elimination of 
u leads to a simple equation for p, its two independent solutions being exponential functions. 
Then substituting p in (jlOj) gives immediately the expression for u. Finally, satisfaction of the 
conditions u = 0, on the two boundaries confining the flow leads to the instability criterion. 
The same procedure is applied here to the cylindrical case. After some calculations, the first 
step gives the equation satisfied by p 
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d-r Lo9 [n) 



p = 



(11) 



In the general case, the solutions of (jlip cannot be expressed in terms of simple analytical 
functions. However, for circular Couette flow with angular velocity given by (JT]) the Rayleigh 
discriminant is proportional to fj and the derivative of the logarithmic term in (jlip vanishes. 
This remarkable feature occurs for angular velocity profiles of the type fi(r) = Cr x + D, 
where C and D are constant coefficients. This allows considerable simplification of Eq. (|11|) 
that, provided j^O, becomes 

'1 n' s 



D z p + 



n 



Dp — m Qp = 0. 



(12) 



It should be noticed that the term in the left-hand-side of Eq. (|12p is reminiscent of the 
potential vorticity describing the content of a Rossby wave as mentioned in [TO]. In the small 
gap limit, rj — > 1, with r/R\ = 1 + (1 — rj)(x/rj) and f2 constant, the term in front of Dp 
disappears from (|12[) thus recovering the solution p = exp(±mv &x) found in [TO] for the 
stratified rotating plane Couette flow when the Rayleigh criterion for stability, $ > 0, is 
satisfied. Then, introducing p = p/vQ, it satisfies 

A / 3A 



£>oi> + 



Or 4 



nr 2 



Am 2 BSl 



p = 



(13) 



with the differential operator C$ = D*D accounting for curvature effects. In Eq. (|13p use has 
been made of the identity m 2 $ = fh 2 Q, where fh = mfli/N and <& = ABCl. We shall assume 
that fh is large enough so that the term in fh 2 dominates in Eq. (|13p . Our analysis differs 
from previous studies that linearized the angular velocity : O(r) = S~i(ro) + (r — vq)£1 + • • • 
around a mean radius r$ as in the thin gap limit. Here, the exact expression of O(r) given in 
dU) is substituted in (|13p that becomes 



A 



C p - 4fh 2 B [ — + B)p = 0. (14) 



When A and B have the same sign, the solutions of (|14p are the modified Bessel functions 
I n (ar) and K n {ar) of general order n real positive with n = 2fh\/AB and a = IfhB. The 
expression for p is 

p = A I n (ar) + B K n (ar) (15) 

where the unknown coefficients Aq and i?o w iH be determined by the satisfaction of the 
boundary conditions on u. Substituting (fi~5|) in ^0]) and using the relation between the 
Bessel functions and their derivatives, one gets u = u/y£l with 

u = A U 1 (r) + B U 2 {r) (16) 

Using vector notations, the two functions U\(f) and ^(r) are considered as components of the 
bidimensional vector U = (Lq,^)- Similarly, introducing the vector *& u = (I u (ar), K u (ar)) 
for v £ (n, n — 1) gives the expression for U 

u = 2 7 $ » +s ((r!!i)*» T H (17) 

In the following, we shall neglect the terms in £1 /CI that appear in factor of a, which are 
small compared to the terms proportional to n or a which behave like fh. 

Satisfaction of the boundary conditions u(R\) = u(R2) = provides the algebraic system 

MA = (18) 

with the vector A = (Aq,Bq) and the elements of the matrix M given by My = Uj(Ri). 
Their expressions are easily deduced from 

U(^)4^T«^Li (19) 

Hi 



where Oi = u> — Cli and Aj = 2fi^ + noi with Jlj = Cl(Ri) for i = 1,2. We have also introduced 
^ n = (In ,Kn ) where the abbreviations In = In(otRi) and if n = K n (aRi) have been 
used. 



2.3 Dispersion relation 

The vanishing of the determinant associated to the algebraic system in Eq. Q18|) leads to the 
dispersion relation 

-5—5-Sn ~ a 2 a 1 a 2 S n -i + aa 2 —Ci 2 - ao 1 —C 2 i = (20) 

ti\ ti 2 ti\ ix 2 

where the following quantities have been introduced : 

S n = I<pKV>-I&>KM and C tj = I® K& + I&K® for (ij) € (1,2). (21) 

To simplify the calculations we shall use the asymptotic expansions at large arguments of 

the modified Bessel functions I and K. At the lowest order, the asymptotic behaviors of 

expressions (|2ip are : 

sinh(q 1 - q 2 ) cosh(ai - q 2 ) , 99 x 

Jn = <-5 n _i ~ , and G12 = 021 ~ , {**) 

^/OL\OL 2 ^/OL\OL 2 

Introducing q = a(R 2 — Ri), the dispersion relation ([20]) reads 

7i - T 2 coth q = (23) 

where 71 = f 2 - a 2 a 1 a 2 J and T 2 = a f a 2 — - a\ -^- J (24) 

are second order polynomials in lj that will be expressed in terms of the rate of shear S given 
below: 

5 = ^— <0, =*► (25) 

^ +1 I ff 2 = © - 1 - 5 

Doing the change of variable w = a) — 1, expressions (j24j) becomes 

ri = X(^-5 2 ) + -i^^ + 5 2 ) + ^-p (26) 

T 2 = £{2[(l-r ? )^ + 5 2 )-(l + r ? )5(l+^)]+n(l-r ? )(^-5 2 )} (27) 

The quantity X involved in expression (|26p for T\ can be expressed in terms of the control 
parameters, \x and 77, as follows 

X = — -a 2 = \rh 2 B^V- (28) 

R1R2 l-r] K J 



2.4 Instability criterion 

It is worth while noticing that the sign of the quantity rj — fi which appears in expression f|28[) 
for X seems determinant for the stability of the system according to what has been observed 
both in experiments [1] and in numerical computations [2J . To check whether our calculations 
support this finding we shall determine the roots of the dispersion relation (|23p . Instability 
could occur if there is a pair of complex conjugate roots, one with negative imaginary part 
corresponding to instability growth, the other to decay. The calculation of the discriminant 
and the determination of its sign is facilitated by introducing the following quantities 

Oi 

t = X - n(l - rj)Q, with Q = — -cothg, (29) 

R\ 

and 

An 
h = — — - 2(1 - rj)Q, fa = 2(1 + rj)SQ. (30) 

The dispersion relation (|23p now reads 

6jH + u(h + fa) - S 2 (t " *0 + *2 + ss-il ~S 2 ) = (31) 

and the associated discriminant is 



A = (ti + fa) 2 - 4fa 



-5 2 (t -fa)+fa + 7? V(l-5 2 ; 

ill/12 



(32) 



Instability could occur if the discriminant is negative. To determine the sign of A two cases 
will be considered according to the sign of fa- Expression (|29|) for fa shows that it is necessarily 
negative when X < 0, which occurs for u > rj. On the opposite side, when X > for /j, < rj, 
fa can takes positive as well as negative values. The change of sign of fa occurs when 

X = n(l - rj)Q (33) 

Replacing X, n and Q by their expressions in terms of n and rj, the above equation becomes 



{1 + rj)(rj - fi) 



cothq = bo, with 60 = 7; ^ — /,., w ». (34) 



Eq. (|34|) can only be satisfied if 60 > 1. Writting that the square of 60 is larger than unity, 
gives the following relation between /x and r/ 

2(1 + V 2 )u 2 - /i[(l + rj 2 ) 2 + 4t? 2 ] + 2(1 + t] 2 )t] 2 > (35) 



The left-hand-side of (|35p is a second order polynomial in \x with coefficients depending on 
r/. It takes a positive value for values of [i ranging outside the interval between the two roots 
/i_ = 2rf /{1 + rj 2 ) and ll + = {l + rj 2 )/2. The root ^+ > rj, introduced when taking the square 
of 60 m (|34p is spurious and will not be considered. Thus, the meaningful condition for the 
satisfaction of 60 > 1 is : 



IX < jX- with jx- = — — ; — 2T (36) 



2rf 
(1+7] 2 ) 
Numerical values of /i_ are reported in Table 1 for different values of rj. The above results 

are summarized as follows. For ll > ll- and bo < 1, we have to < 0, whatever the value of q. 

For fi < /X-, and bo > 1, one can write bo = cothqo and the sign of to depends on the value 

of q. One gets 

to > when q > qo or to < when q < qo- (37) 

Considering separately the two cases, to < and to > 0, we shall derive the conditions for 
instability (A < 0) in the next sections. 

2.4.1 t < 

It was shown in the previous section that to can be negative for any value of n but with 
the additional constraint cothg > 60 when fi < /i_. When to < 0, a necessary condition for 
instability could be derived from the requirement that the term in factor of to in Eq. (|32|) is 
negative 

-5 2 (t -t 1 ) + t 2 + -i-(l-5 2 )<0 (38) 

Replacing to, t\ and ti by their expressions (|29]) - ([30]) and proceeding to some manipulations, 

one gets 

- S 2 X + SQ [2a + nS(l - 77)] + -^-5- [(1 - S 2 ) + nS 2 ] < (39) 

R1R2 L - 1 

where the positive quantity a is expressed in terms of /x and rj 

a = l-S + 7 ] (l+S) = 2il + m) <2 (40) 

r i [X 



Then, expression (|28p for X is substituted in the left-hand-side of Eq. (J39J) that is written as 
the sum of two terms, leading to 

Pi + S 2 P 2 < 0, (41) 



10 



with 

Pi = a 2 S 2 + 2aS^-cothq + —t- (42) 

H\ K1K2 

( n _ 2) 2 a 

Ps = ~ /? p +(l-7?)n— cothg (43) 

It is worth noticing that in the small gap limit (77 — ► 1) the term P2 is irrelevant since the 
terms in n do not exist and S 2 « 1. In this limit, when R\ ~ P2 and a ~ 2, Eq. (|39[) 
reduces to : Pi < 0, where Pi can be factorized, to get 

aS + — - coth I ) f aS + — tanh | ) < (44) 

Pi 2/ V Pi 2/ 

thus recovering the condition for instability, S < 0, found for stratified rotating plane Couette 
flows. The instability occurs for values of q that belong to an interval bounded by the two 
values, q- and q + , deduced from 

tanh — = 7 — , and coth — = 7 — with 7 = 7 r. (45) 

2 ' 2 ' 2 ' 2 ' (I-77) v ' 

For 7 > 1, the instability region is q < g+ while for 7 < 1 the values of q are in the interval 
g_ < q < q + . It is found that 7=1 when fi = 2r/ — 1. For r/ < 0.5, the value of 7 is lower 
than unity, whatever the value of /j. 

When curvature effects are present, the occurrence of instability could be suppressed when 
the term Pi is positive and 5 2 p2 > |Pi|. As it will be intricate to derive a global condition 
for satisfaction of Eq. (|4ip it will be replaced by a stronger constraint that consists to impose 
separately Pi < and P2 < 0. Assuming that (n — 2) 2 ~ n 2 in Eq. (|43p . then P2 will be 
negative for values of q satisfying the condition given below 

61 n /1-^y/ 2 

;i^y where 6l = ^ = H^? 



coth g < y^ where 61 = — = r/ ( ^ ) (46) 



The upper bound for coth q in (|46p has to be larger than unity, this occurs for 

/j^/x with /i = 7?2 + (1 _ 7?) 2 ( 47 ) 

For values of /j lower than /j*, P2 will be negative for values of q such that q > q* with 
cothq* = 6i/(l — 77). We have checked for some representative values of \x and ij that values 
of q inside the interval [g_,g+] are larger than q*. Therefore, the two conditions Pi < and 
P2 < can be satisfied simultaneously. The stability limit fj,* admits asymptotic values. In 

11 



the narrow gap limit (r/ — > 1) it is found that n* — ► 1, which is the stability line for plane 
Couette flows, while in the wide gap limit 11* — > 2t/ 2 (1 + 77). Moreover, the value of the gap 
size for which /j,* = 77 corresponds to 77 = 0.38. Numerical values of ji* are reported in Table 
1 for different values of 77. 

Eq. (|47p provides a necessary condition for instability. To be a sufficient one the positive 
term [t\ + t-^) 1 in expression (|32p for the discriminant should not exceed the negative term. 
The term (t\ + £2) could even vanish if 

coth q = h y ^4 = coth 912 (48) 

(1 -fir]) 



The value coth qi2 is consistent with the values of coth q allowed by Eq. (j46j) if 

(1+f<) <7^ => K* (49) 



(1 - fir]) (1 - 77) 

When 77 > 0.38, the condition, 11 < 77, provides a stronger constraint than the condition 
/i < ]i* coming from Eq. (|47p . 

We have not found the specific value, ij, = (j, s , for which the positive and negative terms 
exactly cancel in the expression for the discriminant written in Eq. (|32p . Thus, the upper 
bound 77 found in Eq. (|49p is a low estimate of the exact value ]i s which is more likely expected 
in the range 77 < jjl s < 11* . 

The above results are stricly valid for values of the rotation rate such that ]X- < 11 < 11* . 
In that case, the constraints on the value of q, which are respectively q £ [g_, gr+] and q > q* , 
can be satisfied simultaneously. When 11 < /j_, an additional constraint has to be satisfied : 
q < qo, as shown in Eq. (|37p . Instability could occur if there is an overlap of the two intervals 
[q_, q + ] and [q* , go] which requires that qo > q-- In that case instability occurs for q 6 [q_, qo\. 

2.4.2 t > 

In the previous section, when investigating the case to < 0, the existence of instability was 
assessed for /i_ < 11 < 77 and q G [q-,q+]- For t/ 2 < /j < /j_, the instability is restricted to 
a narrower range of values of q and its existence was not demonstrated in a systematic way. 
For fj, < IX.., we shall now investigate the case to > that could be more favorable to SRI. 
When to > 0, the terms in the left-hand-side of Eq. (f32j) are rearranged to read 



A = (t 1 + t 2 -2St ) 2 + A 2 (l-S)t t 3 with t 3 =s(-^-Q)-^±p- (50) 

\K\K2 / ti\lX2 
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In that case, the discriminant could be negative if £3 is negative. The change of sign of £3 

occurs when 

coth q =±(^ + ^±p-)=b l + b 2 (51) 

H2 \a a\o\ ) 

The value of coth q in Eq. (|5ip is the sum of two positive contributions. The first contribution 

coincides with b\ given in Eq. (|46p . it is independent of fh and for /j, < [i- it stands in the 

range : 1 < b\ < bo, so that we can write b\ = coth q\. The second contribution 62 ~ a -1 , 

that behaves like fh~ l can be neglected for large values of fh. In that case, £3 will be negative 

for values of q such that 

cothg<6i or q > qi, (52) 

which is consistent with the condition : cothg < bo, allowing for £0 > 0. As soon as the value 
of fh decreases, the contribution 62 increases until it reaches a value such that 61 + 62 = 60 . 
In that case, £3 will be negative for values of q corresponding to coth q < bo which exactly 
coincides with the condition ensuring £0 is positive. Therefore, considering smaller values of fh 
will never give a stronger condition and moreover it will be contradictory with the assumption 
leading to Eq. (fT3|) . 

When £0 > 0, the necessary condition for SRI is of the type 

coth q < coth q max with q < q ma x < Qi, (53) 

meaning that instability might occur for q > q m ax, the value of q ma x depending on the value 
of in. To find a sufficient condition for SRI it will be argued as in the previous section where 
we have looked for the value of coth q corresponding to the vanishing of the positive terms 
in the discriminant. In Eq. (I50p it occurs when £1 + £2 = 25£o, that could be solved only if 
£1 + £2 < 0. Having determined in Eq. (J48]) the value of coth q = coth qi2 for which the sum 
ti +£2 vanishes, the sum will be negative for cothq > coth 512 or correspondingly q < qw As 
coth (712 > cothgi, and consequently qi2 < qi, the simultaneous satisfaction of q > q ma x and 
q < qu is only possible when qmax belongs to the interval qo < q m ax < Q12 which implies to 
consider values of fh not too large. In that case, instability will occurs for q ma x < Q < Qi2- 
Hence, whatever the sign of £0, the values of q leading to instability are restricted to a limited 
band. 



13 



2.4.3 Comparison with experimental and numerical results 

The experimental results pQ obtained for r\ = 0.8 are in qualitative agreement with a first set 
of numerical results [2] obtained for r] = 0.78. Both studies show that the strato-rotational 
instability occurs for /i < rj, in agreement with the condition found in Eq. (|49p . However, quite 
recently we have had knowledge of a second set of numerical results [3] that do not confirm 
these previous findings. The stability line found in ^3] is better represented by /j, = /j,*, where 
/i* is given explicitly in Eq. (|47p . We have reported in Table 1 some values of fj,* corresponding 
to values of r\ used in experiments and computations. In the narrow gap limit, the asymptotic 
value fj,* = 1 is found to agree with the stability line for stratified plane Couette flows. The 
computations performed in [3] for three values of the gap size and different values of the 
Froude number (0.5 < Fr < 2.2) exhibit stability lines /i = /x s with values of jjl s above the 
line \x = T) for r/ = 0.78 and r] = 0.5, in disagreement with numerical and experimental results 
obtained earlier 0[T] for respectively r\ = 0.78 and r/ = 0.8. For these values of 77, the stability 
lines in [3] are in the range 77 < /x s < /i* , in reasonable agreement with the results in Eq. (|47p . 
In the wide gap limit, the necessary instability condition (|47h takes the asymptotic form 
: /i < 2r/ 2 (l + 7]), that fits with a good accuracy the numerical results found in [3] for 
r] = 0.3. The behavior of the numerical neutral curves (Reynolds number versus /i) is strongly 
dependent on the azimuthal wavenumber [2] and on the Froude number [3]. These features 
are not be reproduced by the present analysis based on Eqs. (i9l)- (fi"0|) which are independent of 
/. Although the Froude number appears in these equations through fh = mFr, the instability 
conditions derived here are independent of fh. 



f] 


"if 


V- 




0.8 


0.64 


0.78 


0.978 


0.78 


0.608 


0.756 


0.971 


0.5 


0.25 


0.4 


0.625 


0.3 


0.09 


0.165 


0.231 



Table 1: Stability limits for centrifugal instability (fi c = rj 2 ) and for SRI (fi = fi*) as functions 
of rj. The line \x = \i- plays a fundamental role to determine the sign of to- 
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3 Conclusion 

We performed an inviscid stability analysis of SRI in a Taylor-Couette system characterized by 
a Froude number, Fr, measuring the relative importance of rotation and stratification. Non 
axisymmetric disturbances were considered with azimuthal wavenumber I satisfying IFr << 1 

and / < Imax- 

Finite gap effects were taken into account more appropriately than in previous inviscid 
approaches. Although assumptions have been made, they never concerned the angular velocity 
profile which is kept equal to Q = A/r 2 -\-B. Thus, gap size effects manifest themselves through 
the quantities A and B which depend on \x and 77, the control parameters of the system. We 
derived a necessary instability condition, fj, < /j*, that fits with a good accuracy the recent 
numerical results of Ref. [3J. A stronger condition, /j < 77, found for rj > 0.38, better fits 
with earlier experimental and numerical results [H [2]. For the small gap value rj = 0.3, the 
numerical results [2j [3] exhibit a stability line /j = fi s where fj, s is in good agreement with 
the asymptotic value //* — ► 2r/ 2 (l +77) found for 77 — ► 0. Unfortunately, in the wide gap limit, 
experimental results are still lacking for comparison. 

The angular velocity profile of circular Couette flow is peculiar since it allows an analytical 
resolution for the pressure perturbations in terms of Bessel functions. This is an essential 
step in the above derivation of the stability criterion for SRI in incompressible fluid. A slight 
change in O(r) can lead to completely different stability results. A flow with constant angular 
momentum (f2 ~ r~ 2 ) obtained when B = in ([1]), was considered in a thin cylindrical shell 
|17j . This type of flow is generally assumed centrifugally stable, its Rayleigh discriminant 
being equal to zero. When $ = 0, Eqs. ([9]) and (fTUj) are decoupled and the above analysis 
for SRI cannot be applied. In that case, the existence of nonaxisymmetric unstable modes 
was proved for unstratified flow in a compressible fluid [17] with an equation of state of the 
type p ~ p 7 . It will be interesting in future work to extend the present analytical approach 
to compressible fluids and to other angular velocity profiles. A first step in this direction 
was achieved in a recent theoretical approach [18] based on shallow-water approximation for 
annular sections of Keplerian disks. 

Acknowledgements 
The author is grateful to the Referees for their pertinent comments and valuable suggestions. 
Thanks to P. Le Gal who brought Reference [3] to the author's attention. 

15 



A Appendix: Limit IFr << 1 



When there is no other restriction on the value of I than IFr << 1, the quantity H = I 2 a 2 — <3? 
cannot be simplified. In that case, the governing equations are 



Z _ 

a D*u H u 

r 

Hu — 1i— p 
r 



i g p 
—i a Dp 



(A.l) 
(A.2) 



that depend on the azimuthal wavenumber I through H. Elimination of u gives the equation 
for p 



a 



^p+II-JTjVp 



n d 



H 



&{fh < H --,)+2--Log\^ 



dr 



tt 



p = 



(A.3) 



which has a structure analogous to Eq. (jlip . the main difference is that H appears instead 
of <J>. The derivative of the logarythmic term is given by 



H_ 



n 



-iW (2fl + ^ 



HQ 



Provided a / 0, Eq. (|A.3p becomes 



d2 p + (1-1t) d p + 



(fh 2 H-%)-2l 2 %-{2£l + a) 



p = 



(A.4) 



The above equation can be simplified by assuming m » I or m » IFr^ 1 . For I = 1 and 
Fr = 0.5 this leads to m » 2. As the critical value of the axial wavenumber is not mentioned 
in [21 [3], it cannot be checked if the assumption is satisfied. After introducing p = H x ' 2 p, one 
gets the following equation for p 



C Q p + fh 2 (l 2 a 2 -$)p = 



(A.5) 



Eq. (JA.5P is a generalization of Eq. (|14p that takes into account the value of the azimuthal 
wavenumber L its resolution is left for future work. 
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